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Abstract 

We use recent data on K~^ —> 7r^e'’'e“, together with known values for the pion form factor, 
to derive the kaon electromagnetic form factor for 0 < < 0.125 (GeV/c)^. The results are 

then compared with predictions of the Linear a Model, a quark-triangle model and Vector 
Meson Dominance. The first two models describe the data at least qualitatively, but the 
simple Vector Meson Dominance picture gives a detailed quantitative fit to the experimental 
results. 

PACS numbers: 13.20.Eb, 13.25.Es and 13.40.Gp 


1. Introduction 

The pion electromagnetic form factor, is well stndied experimentally. Many measnre- 

ments, for both positive and negative have been reported in the literatnre [1]. Additional 
information comes from the pion charge radius, which is related to the slope of the form factor at 
= 0 (see [2] for a summary of measurements of the pion charge radius). By contrast, much less 
is known for the kaon form factor, There are some measurements [1] for negative but 

no data exist for > 0. 

Information on Fk can be deduced indirectly from experimental data on the decay —> 

7r+e+e“, such as that provided by the recent high-statistics Brookhaven experiment E865 [3]. The 
amplitude for this decay was measured for g^ up to 0.125 (GeV/c)^, the maximum allowed by the 
kinematics of this kaon decay. However, this amplitude does not give Fk directly, but rather the 
difference Fx^q^) — iU(g^). Since F^{q‘^) is relatively well known, then, this decay is a source of 
information on Fk. 

The Linear Sigma Model (LaM) of Gell-Mann and Levy [4, 5, 6] has proved to be fruitful in 
describing low-energy properties of mesons. In a recent theoretical paper [5], the LaM, and some 
related models, were used to study meson electromagnetic form factors. The authors restricted 
their comparison of iU and Fk with experiment to a check of the slopes at g^ = 0, using the 
measured pion and kaon charge radii. Here, we extend the comparison of these models with 
experiment, and we examine their success in reproducing Fi^(g2) for 0 < g2 < 0.125 (GeV/c)^ 

In section 2, we extract the kaon form factor from the experimental data for K~^ 7r+e^e“. In 

section 3, we derive both the pion and kaon form factors from the LcrM. In section 4, we introduce 
the quark-triangle (QT) model and compare the QT predictions with experiment and in section 
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5, we give the Vector-Meson-Dominance (VMD) model predictions. In section 6, we examine 
the predictions of these models for the slopes of the form factors at = 0. Finally, section 7 
summarises our results. 


2. The decay 7r~^e~^e and the kaon form factor 


The decay K~^ 7i~^e~^e~ has been studied theoretically for many years. Already in 1985 it 

became clear [7] that the process is dominated by the “long-distance” (LD) terms, in which a 
virtual photon is radiated by either the pion or the kaon. However, it was not until the detailed 
data of experiment E865 [3] became available that a convincing description of both the scale and 
the dependence of the amplitude was found [8]. 

Burkhardt et al. [8] considered four contributions to the amplitude, depicted in hgure 1. The 
LD terms are those in figure 1(a) and (b). These two graphs are related to the pion and kaon form 
factors as shown below. Figure 1(c) represents all short-distance (SD) terms. These were already 
known in [7] to be small, and subsequent work [9] has shown them to be smaller than previously 
believed. Therefore here, as in ref. [8], we neglect the SD contribution from figure 1(c). Figure 
1(d) is a “pion loop” term, hrst discussed by Ecker et al. [10]. Its contribution is small, but it 
gives a characteristic shape to the dependence of the amplitude. As in [8], we take this term 
directly from [10]. Burkhardt et al. [8] gives a more detailed discussion of the contributions to 
the amplitude. 

The pion and kaon form factors enter via the graphs of hgure 1(a) and (b), which give the LD 
amplitude [8] 
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In the numerical calculations below, we take the value of the weak matrix element (tt’*' | \ K^) 

from [11]: 


I (7r+ I Hw I 77+) 1= (3.59 ± 0.05) x 10“® GeVl 
Adding the pion loop amplitude, hgure 1(d), from Ecker et al. [10], we obtain 
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Figure 1: Graphs for n~^e~^e~. (a) and (b) are long-distance graphs, (c) is a short-distance 

graph and (d) is the pion loop term. In each graph, the blob denotes the weak (strangeness¬ 
changing) vertex. The (off-shell) photon converts to 

To apply equation (3), we need experimental values for A(g^) and F,r(Q'^)- For the former, we 
use data from Brookhaven E865 [3]. Their amplitude f{q^) is related to our A{q^) by 

I A(q^) 1= (4) 

where Gp is the Fermi constant and a is the hne structure constant. 

The experimental values of Ft^ from [1] are not in general at precisely the required values of 
q^. However, the data, which are plotted in hgure 2, are well described by the VMD model using 
a rho-meson pole: 



for the region of q^ of interest. This is shown by the solid line in hgure 2. Therefore we use 
equation (5) to calculate the required values of We emphasise that the values we calculate 
this way are essentially experimental; the VMD prediction is used basically as an interpolating 
function between the measured experimental points. 
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Figure 2: Pion electromagnetic form factors squared. The points are the experimental data. The 
curves are for VMD (solid line), LaM (dashed line) and QT (dotted line). The dashed and dotted 
lines are practically indistinguishable in the hgure. 


The relative sign of Ald and A-^ioop was already established in [3] and [8]. To derive Fk from 
equation (3), we need also to determine the sign of Fk — TV- To do so, we observe that in all 
models discussed below, as well as in the available data for < 0, Fk differs less from unity than 
does F^, i.e. 

I ~ 1 I < I -^TT — 1 I • 

This dehnes the required sign, giving 


Fr — Ft, 


q^{m\+ 


mz 


e2 I (7r+ I Hw I K+) 


A{q^) - AT,ioop{q^) 


( 6 ) 


The extracted values of | Fk P are listed in table 1 and plotted in hgure 3 which also shows 
previous data for < 0. The errors in our values of Fk arise from experimental errors in the 
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Table 1: Experimental values for the kaon form factor from the present analysis. 


q^ 

{GeV/cf 

1 FK{q^) V 

q^ 

{GeY/cf 

1 FK{q^) V 

0.0244 

1.071 ±0.083 

0.0744 

1.242 ±0.082 

0.0294 

1.086 ±0.083 

0.0794 

1.257 ±0.082 

0.0344 

1.102 ±0.083 

0.0844 

1.280 ±0.081 

0.0394 

1.118 ±0.083 

0.0894 

1.300 ±0.081 

0.0444 

1.134 ±0.083 

0.0944 

1.316 ±0.081 

0.0494 

1.150 ±0.082 

0.0994 

1.344 ±0.081 

0.0544 

1.169 ±0.082 

0.1044 

1.364 ±0.081 

0.0594 

1.187 ±0.082 

0.1094 

1.396 ±0.081 

0.0644 

1.204 ±0.082 

0.1144 

1.413 ±0.081 

0.0694 

1.222 ±0.082 

0.1194 

1.431 ±0.082 


amplitude A{q‘^) for K~^ 7i^e~^e~ and the error in the weak matrix element (tt^ | H\y \ K^), but 

predominantly from the errors on the experimental values of The VMD model, equation (5), 
gives a reasonably unambiguous value for However, its use is only justihed to the extent that 
it agrees with the experimental points in hgure 2. The weighted RMS deviation of the positive-g^ 
experimental points in hgure 2 from the VMD prediction is 0.0837, and we take this as the error 
in I |T 


3. The Linear a Model 

The quark-level Linear a Model (LaM) has been discussed in several papers (e.g. [4, 5, 6, 12, 
13]). Scadron et a/.[5] gives the LaM expressions for meson form factors in the chiral limit (CL), 

F{q^)LaM = -4i/iVe/(167r^) dx J d^p + x{l - x)q^]-^. (7) 

Here, Nc, the number of colours, is related to the meson-quark coupling, g, by [5, 12] 

g = ‘I'Kj Y^/V. (8) 

Thus equation (7) becomes 

F{q^)LaM = -^ [ dx [ [/ -m^ + x{l - x)q^]-‘^. (9) 

TT^ JO J 

In equations (7) and (9), setting m equal to the mean up-down constituent quark mass, m = 
{ttIu + 'md)/2 gives the pion form factor, iV- For the kaon form factor, m is set to the mean of m 
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Figure 3: Kaon electromagnetic form factors squared. The solid points are the experimental data 
from the present analysis and the circles show the previously existing data. The curves are for 
VMD (solid line), LaM (dashed line) and QT (dotted line). 

and the s quark mass, rUus = + m)/2. In the present paper we use values in the chiral limit 

[14]: 


= 325.7 MeV, = 469.0 MeV, = 397.4 MeV. (10) 

To avoid the divergence in equation (9), we use a Taylor-series expansion as follows. Differen¬ 
tiating equation (9) gives 


2i f d*p 

dq^ TT^ Jo J — mJ + x{l — x)q^]^ 

The second integral in equation (11) can be written 


( 11 ) 
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d^p 

[p2 - l/jS 


( 12 ) 


I = 


where V = m? — x{l — x)q^. Applying the Wick rotation, k‘^ = — d'^p = in‘^k‘^dk‘^, this becomes 


I = 

which can be integrated directly to give 


de 

'-F - l^]3 


and hence, from equation (11), 


I = 


in 


2 


21^ 


(13) 


(14) 


= /■ (15) 

dq^ Jo V 

As a check, the integral in equation (12) can be evaluated from the expression given by Scadron 
[16], and also by integration of the Taylor series for the charge radius given in [8]. 

At = 0, y = gives 


dF{q 


2\CL n 
ILaM 


dq"^ 


J<?2=0 



x)dx 


1 

Qm? 


(16) 


Further derivatives follow from differentiating equation (15). In writing a Taylor series for F{q^), 
the first term, F(0)£^jy^, is given directly by the normalisation requirement that F[q^) = 1 at 
= 0. The series, therefore, is 


F{q 


2\CL 

ILo-M 


= 1 + 


1 M 

6 \m 
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60 Vm 


4 1 

+ 


420 


(-) 


6 1 
+ 


2520 Vm 


Q 


+ 




13860 \m 


10 


+ .. 


(17) 


As in equation (9), this series gives F^^i^q^) with m = rh and Fx^q^) with m = nius- 

The form factors FT,{q‘^)2^j^ and FK{q^)^^M predicted by equation (17) are plotted in hgures 
2 and 3, together with the data from the present analysis and from [1], For both FT,{q'^)2^j^ and 
FKid^L^, the LaM shows the same qualitative trend as the data, but falls somewhat below the 
experimental points. However, equation (7), from [5], is derived using quark loops only. Several 
papers (e.g. [17, 18]) have studied the addition of meson loops to these. Although they do not 
affect the result at g^ = 0, they may be important away from this point, which may be the source 
of the discrepancy. 
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Figure 4: Quark-triangle graphs. 

4. Quark triangle graphs 

Scadron et al [5] describe another approach to meson form factors, via the quark-triangle 
graphs of hgure 4. From [5], the pion and kaon form factors are given by 


F^{q‘^)QT = -^ig^Nc + , 


FK{q‘^)QT = -^ig'^Nc + , 


where 


I{q^,ml,ml,M^) = 


F F -I-2(M^ — (mi — m2)^)(l — m) 

2{27i)'^ Jo ^ Jv ^ ml — {M"^ + ml — mi)u + + {v"^ — u^)q‘^/A 


+ J dx J d'^p p^ — ml + (M^ + ml — m\)x — 


Here, we replace m„ and m^ by m from equation (10). These expressions can be calculated either 
by direct evaluation of the integral in equation (20) or, as a check, by construction of a Taylor 
series by repeated differentiation of equation (20). 



The results for the pion and kaon form factors are shown in figures 2 and 3. As for the LaM, the 
predictions for both Fk and follow the data qualitatively, but detailed quantitative agreement 
is lacking. As with the LaM, the problem may lie with the neglect of meson loops, which could 
be signihcant when > 0. 

5. Vector meson dominance 

Finally, we examine a simple model based on VMD. As discussed in section 2 above, the pion 
form factor is dominated in this model by the p pole, equation (5), and hgure 2 shows that this 
simple picture agrees well with the data. 

For the kaon form factor, there are contributions from the p, uj and 0 poles: 


Fk{q‘^)vmd — N 


1 9 pee 1 9uJee /I 9<j>ee \ 

2 m^o — 2 V 2 )' 


( 21 ) 


where Ppee = 4.97, g^^e = 17.06 and g^ee = 13.38, derived from the decay widths. The p^K^K~, 
ooK^K~ and (j)K~^K~ SU(3) coefficients are 1/2, 1/2 and l/sqrt2 respectively. The requirement 
that F(0) = 1 gives the normalisation coefficient as N = 0.03682 GeV^. 

The prediction of equation (21) is plotted with the data for Fk in figure 3. As for F^,, the 
VMD gives excellent agreement with data, both for the previously available data for < 0 and 
for the new data derived in the present paper. 


6. Slopes of form factors at g^ = 0. 

In this section we examine two properties that depend on the slopes of the form factors at 
g^ = 0, i.e. the meson charge radii and the amplitude for T:^e^e~ at g^ = 0. Predictions 

for charge radii have been published before [5] but we include them here for completeness and 
because we update some numerical values. 

The charge radius is related to the form factor by [5] 


r = 



\ 


dF(q^ 

dq"^ 


Jg2=0 


( 22 ) 


For the LaM and VMD, the quantity dF{q^)/dq^ is straightforwardly obtained from equations 
(16), (5) and (21). For the QT model, the expression is given by [5] as 


and 





'hcY 

.m. 


( 23 ) 
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Table 2: Meson charge radii and amplitudes for at = 0. 


Model 

r-jT 

(fm) 

r-K 

(fm) 

A(0) 

10-9 GeV-2 

LaM 

0.606 

0.497 

7.56 

QT 

0.606 

0.514 

6.47 

VMD 

0.623 

0.574 

3.69 

Experiment 

0.672 ±0.008 

0.56 ±0.03 

4.00 ±0.18 


(4.) ^ 


5 ^ 3 ^2 4 ^3 22 ^4 2 . , 

-5 + -5^ - -S^ H--(5® + ... . 

6 5 9 63 7 ' 


dTT^m^ 

where 6 = {ms/m) — 1 = 0.44 (see also [19]). 

The predictions are compared with the experimental values from [2] in table 2. 


(24) 


The LD contribution to the amplitude for n^e^e at = 0 is given by equation (1) as 


1 Ai^(O) 1= 

(7r+ 1 Hw 1 K+) 

2 2 



^K+ - ^n+ 



dF^ 


dFj^+ 


dq^ dq^ 


g2=0 


(25) 


Since the pion loop term vanishes at q^ = 0, equation (25) gives the total amplitude directly. For 
the three models considered here, A(0) can be calculated using the expressions for the derivatives 
given above. We take the experimental value from Brookhaven E865 [3], converting their /o to 
our A(0) using equation (4). 

From table 2, we see that the VMD picture again comes closest to a quantitative description 
of the data. As expected, this is consistent with the conclusion from the form factors of hgures 2 
and 3. 


7. Summary 

We have derived new experimental values for Fk for positive q^ from the decay Tr+e+e , 

extending the upper limit of the range over which experimental values of Fxiq^) are known to 
q^ = 0.125 (GeV/c)^. Also, we have calculated both Ft, and Fk for the range of q^ covered by 
these new data and the previous data, and have examined the slopes at q^ = 0. It is apparent 
that for all these comparisons, the LaM and the QT model give at best a qualitative description 
of the data, probably due to the neglect of meson loops. By contrast, the VMD model gives a 
detailed quantitative £t both to the data for F,, and to the new data for Fk- 
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